Lattice QCD study on if/* (892) meson decay width 
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We deliver an exploratory lattice QCD examination of the if* (892) meson decay width with the 
help of the p-wave scattering phase Si of pion-kaon (ttK) system in the isospin 1=1/2 channel, which 
is extracted by the modified Rummukainen-Gottlieb formula for two-particle system with arbitrary 
masses, and it clearly reveals the entity of a resonance at a mass around if* (892) meson mass. The 
effective range formula is applied to describe the energy dependence of scattering phase and we 
obtain the effective if* => ttK coupling constant as gK*^K — 6.38(78), and subsequently achieve 
the decay width as Tk* = 64.9 ± 8.0 MeV, which is in reasonable accordance with the experiment. 
Our lattice investigations are conducted on a 20 3 x 48 MILC full QCD gauge configuration at 
(m n + mfc) /mK* ~ 0.739 and the lattice spacing a « 0.15 fm. 

PACS numbers: 12.38. Gc, 11.15.Ha 



I. INTRODUCTION 

It is well-known that the vector if* (892) meson is a 
resonance. In 2012, the Particle Data Group (PDG) 
listed the if* (892) meson I(J P ) = 1(1"), with a mass 
891.66±0.26 MeV and a narrow width 50.8±0.9 MeV Q. 
Some recent experimental analyses have precisely 

measured its resonance parameters. Moreover, a few the- 
oretical efforts have been taken to calculate its hadronic 
coupling constant 0-0 • Since the if* (892) meson is a 
low-lying vector meson with strangeness, a study of its 
decay width is definitely a straightforward investigation 
on the three-flavor structure of the low-energy hadronic 
interactions, hence, it is very helpful for us to acquaint 
with the dynamical traits of the hadronic interactions 
with QCD. 

At present, lattice QCD is the most feasible approach 
from first principles to extract the resonance parameters 
of if* (892) meson nonperturbativcly. The principal de- 
cay channel (with a branching rate of 99.9%) of if* (892) 
meson is to one pion plus one kaon in the p-wave [l[, 
which can then be precisely dealt with on the lattice, and 
there is a pioneering lattice QCD study on its coupling 
constant qk'-kK through evaluating appropriate three- 
point correlation function [Toj . Among unstable hadrons, 
the vector p meson is ideal (see reasons in [ll[) for lat- 
tice QCD inve stig ations of a resonance, and it is exten- 
sively studied [RHllI], nevertheless, so far, lattice QCD 
research on the resonance parameters of if* (892) meson 
directly from p-wave scattering phase of 7rif system in 
the 1=1/2 channel has not been reported yet, mainly 
because the rectangular diagram is exceptionally hard to 
rigorously calculate, the statistical error of the numcri- 
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cally computed if* mass is not too small, and most of 
all, the proper finite size formula to describe 7rif sys- 
tem enclosed in a cubic box in the moving frame is not 
completely established yet until recently. 

Motivated by the recent extensions and developments 
of the Rummukainen-Gottlieb formula [l9| to a generic 
two-particle system with arbitrary masses in the moving 
frame (20l-|26| and J. Nebreda and J. Pelaez's brilliant ex- 
positions on if* (892) resonance [13], and also encouraged 
by our previous work on the accurate determination of 
if* mass [28| , the exploratory calculations of the scalar 
meson decay widths [2lL [29l |30| , and the accurate compu- 
tation of the s-wave 7rif scattering length in the 1 = 1/2 
channel [3l|, here we will step out further to probe its 
decay width by way of lattice QCD. 

In the present work, we discuss all the possible com- 
putation scheme for calculating the 7rif scattering phase 
with total zero momentum in the center-of-mass (CM) 
frame, and total non-zero momentum in the moving 
frame (MF), respectively, and obtain the if* (892) decay 
width by calculating p-wave scattering phase of 7rif sys- 
tem in the 1=1/2 channel in the moving frame. The cal- 
culations are launched on a MILC full QCD gauge config- 
uration with the 2+1 flavors of the Asqtad improved stag- 
gered quarks 132}, [3311 . The meson masses quoted from our 
previous work [28[ yielded (m^ + m^/rriK- ~ 0.739, and 
the lattice ensemble parameters determined by the MILC 
collaboration gave the lattice extent L rj 3.0 fm and 
the lattice space inverse 1/a = 1.373 GeV [HI, [33|. The 
Liischer formula [H|-[36[ is, as a usual, applied to the case 
in the center-of-mass frame, and we use a newly estab- 
lished finite size formula, which is the generalization of 
Rummukainen-Gottlieb formula [l9[ to the generic two- 
particle system in the moving frame [20l - [26| . to estimate 
the p-wave 7rif scattering phase in the 1=1/2 channel. 
The simulations conducted at two energies around the 
if* resonance mass enable us to extract the decay with 
of the if* (892) resonance. 
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This article is organized as follows. In SccQll we elab- 
orate on our calculation method. Our concrete lattice 
calculations are provided in Sec. Mil We deliver our re- 
sults in Sec. IIV1 and reach our conclusions and outlooks 
in Sec. El Numerical calculations of the zeta function are 
courteously supplied in the appendix for reference. 



B. Finite- volume methods 

In this paper, we only focus on the irK system with 
the isospin representation of (I, I z ) = (1/2, 1/2) and de- 
liberate on the if* (892) meson decay into one pion plus 
one kaon in the p-wave. 



II. FORMALISM AND METHOD OF 
MEASUREMENT 

A. The relativistic Breit-Wigner formula 

The -KT*(892) resonance possesses quantum numbers 
I{J P ) = ) and principally decays into one pion and 
one kaon in the p-wave with a branching rate of 99.9% [l[. 
For an elastic ttK scattering in the resonance region, the 
relativistic Breit-Wigner formula (RBWF) for the p-wave 
scattering phase ft can be written as [l[ 



tan ft = 



\/sTr{s) 
M%-s ■ 



s = E, 



(1) 



where Mr is the resonance position, Y r is decay width, 
Ecm is the center-of-mass energy, and s is the Mandel- 
stam variable. The Tr(s) can be expressed by way of the 
effective K* — > irK coupling constant qk-^k as [271 ] . 



r*(*) 
p 



9k*-kK P 3 
6tt s 



(2) 



2^ 



y/[s - {m-n - m K ) 2 } [s - (rrhr + m K ) 2 ], 



Checking equations (JXJ) and © , a representation of the p- 
wave scattering phase as a function of the invariant mass 
yfs is offered by the effective range formula (ERF), 



tan ft = 



6ir V^(M R - s) ' 



(3) 



which is applicable in the elastic region and suits the 
experimental measurements pretty well. The ERF per- 
mits us a fit or seeking for two unknown quantities: the 
coupling constant qk'^k and the resonance position Mr 
from the p-wave scattering phases. Then the K* decay 
width is computed by 



K* 



Pk* 



=Ml 



9k*itK Pk* 
6tt M% ' 



(4) 



2M R 



Equations (J3)) and ((4]) provide us an approach to derive 
the decay width r k* by studying the dependence of the 
p-wave ttK scattering phase Si on the invariant mass y/s. 
We should stress at this point that we will extensively 
apply the ERF approximation in the present study since 
the RBWF holds perfectly for relatively narrower objects 
and the K* (892) resonance has a pretty narrow decay 
width 50.8 ±0.9 MeV Q. 



1. Center of mass frame 

In the center-of-mass frame, the energy eigenvalues of 
the non-interacting ttK system reads 



K 



+ P, 



where p = |p|, p = (27r/L)n, and n g Z 3 . The energies 
for the n^O are typically larger than the K* resonance 
mass rriK* ■ For example, the lowest energy for the n ^ 
calculated from the previous determinations of m ff , my 
and m,K* IH} is E as 1.12 X mx*, which is self-evidently 
not qualified to study the K*(892) meson decay. Hence, 
we have no choice but to consider the n = case, and 
the energy E = 0.739 x mjf«, which is still not a favorite 
option. 

When considering the interaction between pion and 
kaon, the energy eigenstates of ttK system are displaced 
by the hadronic interaction from E to E, which are cal- 
culated by 



E = s/ml + k 2 + 



k = 



2tt 
1 



'I- 



where the dimcnsionless momentum gel. Solving this 
equation for the scattering momentum k, we have 



k 



1 

2! 



[E - (niir - m K ) 2 ][E - (m^ + m K ) 2 



In this article, we are primarily interested in the energy 
eigenstates of ttK system in the elastic region m^+rriK < 
E < 2(m- K + mx). In the center-of-mass frame these 
energy eigenstates transform as a vector (to be specific, 
the irreducible representation T = T±) under the cubic 
group Oh ■ The p-wave ttK scattering phase Si is linked 
to the energy E by the Liischer formula |33 - l3a |. 



tanft(fc) = 



r 3/2 



q 



Zoo(i;q 2 )' 

where the zeta function is formally defined by 

1 x - 1 



2 o(s;q 



E 



n 



- q 2 Y 



(5) 



(6) 



The Zoo(s;q 2 ) has a finite value only when Res > 3/2, 
nevertheless it could be analytically continued to s = 1. 
We usually evaluate Z o(s;q 2 ) using the way described 



in Ref. 37[. I notice that there exists an equivalent 
Liischer formula in Ref. |38j |. which is the generalization 
of the Liischer quantization condition to multiple two- 
body channels. Moreover it is easy to calculate and more 
accurate than Liischer formula in the relativistic case. 
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2. Laboratory frame 

To implement the physical kinematics such that the 
energy of irK system is pretty close to K* meson mass, 
we recourse to the laboratory frame [l9| , which is usually 
called the moving frame. We have presented the detailed 
discussions of irK system in the moving frame in Ref . [30| , 
here we just review its essential parts. 

Considering a moving frame with non-zero total mo- 
mentum P = (2w/L)d, d € Z 3 , the energy eigenvalues of 
the free pion and koan are 



Emf = \j ml + pj + ^m 2 K + p z 2 , 

where pi = |pi|, P2 = |p>2 1 , and pi, P2 define the three- 
momenta of 7r and K, respectively, which meet the peri- 
odic boundary condition (PBC), 

2tt 2tt , 

Pi = — n x , p 2 = — n 2 , ni,n 2 eZ'\ 

and total momentum P satisfies P = pi + p 2 . 
In the center of mass frame, the energy Ecm is 



E, 



CM 



P 



P 



where total center-of-mass momentum disappears, 
namely, p* = |p*|, p* = pj = — p 2 , here we denote the 
center-of-mass momenta with an asterisk (*) [l9j]. We can 
readily verify that the p* are quantized to the values [30| 



2tt 

P = — r • 
V L 



where the set Pd is 

Pd = <; 



d 

2 



E 2 

^CM 



n € Z c 



(7) 

where the boost factor, 7 = l/Vl — v 2 , operates in 
the direction of the velocity v, which is calculated from 
v = P/ Emf 7 an( l for the notational compactness we have 
taken the shorthand notation Oil , 



7P = 7PII +P-L, 



7 *P 



7 ^ll 



P±, 



(8) 



where py and pj_ are the ingredients of p parallel and 
perpendicular to the center-of-mass velocity v, respec- 
tively [H: P|| = (p • v)v/|v| 2 , p ± = p p . 

Using the standard Lorentz transformation, the energy 
Ecm is connected to the Emf through Ecm = j~ 1 Emf, 
or by EIm = E MF - P 2 . 

We are particularly interested in one moving frame: 
pion at rest, kaon with the momentum p = (27r/L)e3 
(namely, d = e^) and K*(892) meson with the momen- 
tum P = p. For our concrete case, we found that its 
invariant mass takes y/s = 0.8788 x niK*, which is sig- 
nificantly closer to K* meson mass rriK* than that in 
the center-of-mass frame. Finally, we find one suitable 
to study K* decay, and here we solely consider this case. 



In the presence of the interaction between pion and 
kaon, the Ecm can be calculated by 



2tt 



E C m = V m l + k2 + V m K + fc2 ' k = ~i q - 

where the dimensional momentum { £ R. Solving this 
equation for the scattering momentum k, we arrive at 



1 ~2 

k = T^V \ e cm - \jrin -m K ) 2 ] [E CM ~ («V + m K ) 2 ] ■ 

(9) 

It is convenient to rewrite equation ((9|) to an elegant form 
for later use as 



1 



'CM 



•K 



E 



CM 



(10) 



which is used to calculate the scattering momentum k 
(including its statistical error), and investigate the lattice 
discretization effect. 

The energy eigenstates of ttK system for our moving 
frame transform under the tetragonal group C± v [2l| . 
Only the irreducible representations A\ and E are as- 
sociated with the p-wave ttK scattering states in a torus. 
We only compute the energies related with the A\ sector 
in the present study. The hadronic interaction displaces 
the energy eigenstate of irK system from E to E, and 
the energy E is linked to the p-wave nK scattering phase 
#i with the help of newly established finite size formula 
in the moving frame for t he g eneric two-particle system 
with arbitrary masses [ITJ, |2l[ |23| , 



77T 3 / 2 <7 



(11) 



where we does not consider the higher scattering phase 
shifts Si(l = 2, 3, 4, • • • ) [221 ]. and the modified zeta func- 
tions are formally defined by 



rGP d 

E 

rGPd 



(|r| 2 -9 2 ) s ' 
r 2 Y 20 (n r ) 
(r 2 — q 2 ) s ' 



(12) 



where f2 r represents the solid angle parameters (6, (f>) 
of r in spherical coordinates and the Y/ m are the stan- 
dard spherical harmonic functions, and the set Pd is 
denoted in Eq. ([7]). The scattering momentum k is 
calculated from the invariant mass ^/s through t/s = 

We have discussed the numer- 



Vk 2 



Vk 2 



K ' 



ical calculation method of the 2^(1; q 2 ) in Appendix 
A of Ref. [3(1, and we will give the numerical calcula- 
tion method of the Zf Q (l;q 2 ) in Appendix [Al although 
there are some general calculations of the zeta function 
Z L^1 2 ) inRefs. HlHl III . 
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Correlation matrix 



To compute two energy eigenvalues, i.e., E n (n = 1,2), 
we constitute a 2 x 2 correlation function matrix: 



C(t) = 



(OIO^ (t)O wJf (0)|0) (0\Ol K (t)O K . (0)|0> 
(0\O^ (t)CW(0)|0) (0\O K , (t)<D K . (0)|0) 



(13) 

where Ok * (t) is an interpolating operator for the vec- 
tor if* (892) meson with the specified momentum p = 
(2w/L)e 3 and the polarization vector parallel to p; 
Officii) is an interpolating operator for the irK system 
with the given momentum p = i2ir/L)e 3 . These inter- 
polating operators employed in the present work are ac- 
tually identical to those in our previous studies [H, H3 > 
nevertheless, to make this article self-contained, all the 
fundamental definitions will be provided in the following 
as well. 



1, ttK sector 

Here we take advantage of the original definitions 
and notations [391443T ] to examine the necessary formu- 
lae for the lattice QCD calculation of the p-wave scat- 
tering phase of irK system enclosed in a torus at the 
I = 1/2 channel. Let us learn the elastic scattering 
of a Nambu-Goldstone pion with zero momentum and 
a Nambu-Goldstonc kaon with the momentum p in the 
Kogut-Susskind (KS) staggered fcrmion formalism. Us- 
ing the interpolating operators 0„(xi), 0^(13) for pions 
at points xi,xa, and Ok(x2), Ok (xi) for kaons at points 
X2,a;4, res pec tively, the -kK four-point functions are ex- 
pressed as|43| 

CffK ix4 ) = (0^x4)0^x3)0^2)0^)). 

where the pion and kaon interpolating field operators are 
denoted by [43] 

O w +ix,t) = -d{x,t)j 5 u{x,t), 

O n o{x,t) = -r=[u{x,t)^ 5 uix,t) - d(x,i)7 5 d(x,i)], 

K o(x,t) = s(x,t)7 5 d(x,t), 
K +(x.,t) = s(x,t)7 5 u(x,t). 

After carrying out the summation over the spatial coordi- 
nates Xi, X2, X3 and X4, we gain wK four-point function 
with the momentum p as [30( 



cw(p;*4,i3,*2,*i) = < 



ip-(x 4 — x 2 ) 



Xi X 2 X 3 X4 

XCff K (x4,X 3 ,X2,Xi), 

where Xi = (xi,ti), x 2 = (x 2 ,t 2 ), x 3 = (x 3 ,£ 3 ), and 
X4 = (x4,^4). To refrain the color Fierz rearrangement 
of the quark lines [4l[, we choose t% 7^ t 2 7^ £3 7^ £4 and 
set t\ = 0, t% = l,t3 = t, and £4 = t + 1, respectively, 



here t represents the time difference. We build the 7r.fr 
interpolating operator in the 7 = 1/2 channel as [43| 



= ^{V2n + it)K (p,t + 



1 



(14) 



where p is total momentum of ttK system. The ttK 
operator has the isospin representation with (J,i») = 
(1/2, 1/2). _ 

Considering that u and d quarks have the equal mass, 
topologically only three quark line diagrams still con- 
tribute to irK scattering amplitudes [IH . These diagrams 
are elucidated in Fig. [TJ and we usually label them as di- 
rect (D), crossed (C) and rectangular (R) diagrams, re- 
spectively [10, 1 . The direct and crossed diagrams can 
be readily computed [4(| 5]J by means of only two fixed 
wall sources placed at the time slices t\ and t 2 , which en- 
ables a relatively chea p la ttice calculation of the 7 = 3/2 
ttK scattering length [44l lilj . Nevertheless, the rectan- 
gular diagram (R) needs extra quark propagator linking 
the time slices £3 and £4, which make the strict evaluation 
of this diagram extraordinarily expensive. 



*4- 
*3- 



(a) D 



(b)C 



(c)R 



FIG. 1: (color online). Quark-link diagrams contributing 
to ttK four-point functions. Short black bars stand for the 
wall sources. Open circles are sinks for local pion or kaon 
operators. The black lines represent for the u/d quark lines, 
and red lines stand for the strange quark lines. 

Sasaki et al. handled this puzzle through the technique 
with a fixed kaon sink operator to save the computational 
resources (46[. Lang et al. recently solved this prob- 
lem by the use of Laplacian-Heavyside smeared quarks 



47] . In our previous 



within the distillation method 

works [2l|, we settled the question using the mov- 
ing wall sources witho ut g auge fixing introduced first by 
Kuramashi et al. [4(| l4l| . More specifically, we calcu- 
lated these diagrams by computing each T quark prop- 
agators corresponding to the wall source at all the time 



1 In Ref. 1431 . they are denoted as A, H, and X, respectively. 

2 It is well-known that the rectangular diagram (or backtracking 
contractions, box diagram [4^1 1 is most challenging and impor- 
tant for the 1=1/2 channel. And obtaining the reliable signal 
of it is vital to our final result. We observe that the signals of 
the rectangular diagram in Ref. [47ll are at reasonable levels. 
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slice [H [H EE E3, 

Y^Dn'^Gtin") = X;*n',(*,t), < t < T- 1, 

Tl" X 

where -D is the Dirac quark matrix, and the subscript t in 
the quark propagator G indicates the position of the wall 
source in the temporal direction [2l], [30, Efl EH . The as- 
sociations of the quark propagators Gt(n) exploiting for 



the ttK four-point correlation functions are schematically 
illustrated in Fig. [1] 0, M E3, S3 ■ As we practiced in 
Ref. for the study of k decay width, for the non-zero 
momentum p, we take an up quark source with 1, and 
a strange quark source with the e lp ' x on each lattice site 
for the pion and kaon creation operator, respectively [30j |. 
By means of the quark propagators G, we can represent 
D, C, and R diagrams as [3(| 



C^(p;t4,i3,*2,*i) = ^^e^(Tr[Gl i (x3,t3)G tl (x3,t 3 )]Tr[Gl 2 (x4,t4)G t2 (x 4 ,t 4 )]), 

C-nK (P! * 3 , h, tl) 
Ckk(Pi * 3 , h, ti) 



EE' 



(Tr^ (x 3 , t 3 )G ta (x 3 , < 3 )GJ 2 (x 4 , U)G tl (x 4 , i 4 )]> 



]T ]T e ipX2 (Tr^ (x 2 , t 2 )G t4 (x 2 , t 2 )Gj 4 (x 3 , i 3 )G tl (x 3 , i 3 )] 

X 2 X 3 



4 



r 



(15) 



where the traces are taken over color, and the hermiticity 
natures of the quark propagator G have been applied to 
remove the 7 5 factors [l^JMJS, S3 . 



In practice, we compute the ratios as well 



As discussed in Refs. 4(1" 41|~the rectangular diagram 
R produce the gauge- variant noise, and we usually reduce 
it by conducting the gauge field avera ge w ithout gauge 
fixing as we practiced in Refs. [2l|, H3, HE EH- All the 
three quark line diagrams in Fig. [T] are needed to compute 
the p-wave irK scattering phase in the 1 = 1/2 channel. 
In the isospin limit, the ttK four-point function in the 
1 = 1/2 channel can be described in terms of only three 
quark line diagrams [43j , 



= d+ 1 -n s c-In s r, 



(16) 



where the interpolating field operator O^k denoted in 
Eq. (TH| generates a irK state with the total isospin 1/2 
and momentum p, and Nf is the staggcrcd-flavor factor, 
which is plugged in to address for the flavor degrees of 
freedom of the KS staggered fermion [39| . We should bear 
firmly in mind that if we carry out the appropriate root of 
the staggered fermion determinant 3 , in the continuum 
limit, the same number of the flavors flow around the 
internal quark loops as in QCD [39|. Therefore, at the 
level of these quark line diagrams (namel y, D , G, and R), 
all contributions are exactly as in QCD [39(. 



iZ*(t) 



C£ f (p;0,l > t,f+1) 



X = D, G, and R, 



C n (0;0,t)C K (p;l,t+l) 

(17) 

where CV(0;0,£) and Cjr(p;l,i + 1) are the it and K 
correlators with the momentum and p, respectively. 

We should bear in memory that the dedications of 
non-Nambu-Goldstone pion and kaon in the intermediate 
states are exponentially reduced for the large time owing 
to their heavier masses compared with these of Nambu- 
Goldstone pion and kaon [HI, H3, 139M41) . Thus, we can 
ignore this systematic error due to other ttK tastes. 



2. if* (892) sector 

In principle, we can calculate the propagators for two 
local vector K* meson, 7.; ® 7^ (VT) and 707* ® 7o7i 
(PV) [(30|, HH . However, here we simply quote the results 
for local VT K* meson since it delivers quite stable re- 
sults in the analysis of the mass spectrum. Moreover, the 
numerical evaluation of K* — > irK three-point function 
is much eased if we adopt local VT K* operator. Thus, 
we use an interpolation operator with the isospin 1=1/2 
and J p = 1~ at the source and sink [28j . namely, 

0(x) = ^ u a(x)ji <8> -ftS a (x), 



3 There are some strong evidences to demonstrate that the con- 
tribution from a single Dirac fermion can be nicely restored by 
carrying out the fourth root of the fermion determinant, see more 
details in Ref [50ll . In the present work, we suppose that the 
fourth-root procedure reproduces the proper continuum limit of 
QCD, and the lattice results of this work rely on this hypothesis. 
Please consult Refs. [5lH59l | for the recent investigations about 
the fourth-root recipe. 



If imposing the Dirichlct boundary condition in the temporal 
direction, we can easily extract the energy shift &E from ratios 



R x [4fj,|4l|]. On the other hand, we can readily check that when 
t <C T/2, even we set the PBC in the temporal direction, we still 
can roughly estimate 8E from these ratios. 
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where a is the color index. The time slice correlator for 
the K* meson in the momentum p state is computed by 

C K ,(p,t) = ^^e ip ' x (u 6 (x,t)7i®7 i s 6 (x,t) 

x a ,b 

xs o (0 J 0)7i®7i«2(0,0)> J 

where 0, x are the spatial points of the K* state at source 
and sink, respectively. 

For the staggered quarks, the meson correlators have 
the general single-particle representation, 

C{t) = Aie- mit + ^ A'ii-lfe-™^ + (t -> N t - t), 

i i 

where the oscillating terms correspond to a meson with 
the opposite parity. For K* meson correlator, we take 
only one mass with each parity, and the oscillating par- 
ity partner is the p-wave meson with the J p = 1 + . The 
K\ meson is with J p = 1 + , so it is the candidate of 
the oscillating parity partner of the vector K* meson. 
However, these states with J p = 1 + can just as well be 
multihadron states 5 . With staggered fermions, the mul- 
tihadron possibilities include the various taste combina- 
tions. So we can not identify its parity partner with the 
Ki, see more discussions in Ref. [28|. Thus, the if* (892) 
correlator was fit to 

C K ,{t) = b K ,e~ m ^ + &x- 1 (-l)*e- M *i* + {t->N t - t), 

(18) 

where bx 1 and bx* are two overlap factors. 

3. Off-diagonal sector 

A calculation of the generic three-point function are 
briefly discussed in Ref. [Io|. To rigorously evaluate it 
we must compute a spatial volume number of propaga- 
tors, namely (16 3 for our case). To avoid the ap- 
parent intractability of the exactly computing this prob- 
lem, S. Golttlicb et al. introduced the "exponential" 
method, which calculate a two-point function with the 
presence of a source, and then differentiates with the 
source strength to achieve the corresponding three-point 
functions 0,11]. To investigate vector meson decay into 
pseudoscalars from quenched lattice QCD Hil. L oft and 
DeGrand adopted "two-stage" technique [62|, [63| , which 
takes approximately twice as compared with the calcu- 
lation of the mass spectra [Ioj |. Later, when stud yin g 
the resonance parameter of the vector p meson |10l4l8|. 
people chiefly employ a stochastic method l(i ll lili or its 
variants to evaluate three-point correlation function. 

Motivated by the precisely evaluate the tttt four-point 
correlation functions by Kuramashi et al [H, |4l[ with the 
moving wall source technique [3lj . analogously, we have 



successfully extended this technique to evaluate three- 
point correlation function, and obtained pretty good sig- 
nals for the three-point functions of the rnr — > a [2!| and 
ttK — > k (3pj . In this work we will continue to use this 
technique to evaluate the irK — > K* three-point correla- 
tion function. 

To prevent the tangled color Fierz transformation of 
the quark lines [4l|, we should choose t\ ^ t<i ^ t^. 
In practice, we pick t\ = 0,^2 = 1, and t% = t for 
the ttK — > K* three-point correlation function, and opt 
t\ = 0,< 2 = t, and t3 = t + 1 for the K* -> irK three- 
point correlation function. The quark line diagrams cor- 
responding to the K* — > ttK and irK — > K* three-point 
functions are schematically illustrated in Fig. [UJa) and 
Fig. Hl^b) , respectively. 




FIG. 2: (color online). Quark-link diagrams contributing to 
ttK — > K* and K* — > nK three-point correlation functions. 
Short black bars stand for the wall sources. The black lines 
represent for the u/d quark lines, and the red lines stand for 
the strange quark lines, (a) Quark contractions of ttK — > K* , 
where open circle is the sink for local K* operator, (b) Quark 
contractions of K* — > nK, where open circle is the sink for 
local pion operator. 

The irK — > K* three-point function can be easily eval- 
uated with only two fixed wall sources [2l|, [3(| . Neverthe- 
less, the computation of K* irK three-point function 
is pretty hard, since it requires the extra quark propaga- 
tor linking time slices t 2 and t 3 [2l|, HI] . In practice, we 
employ an up quark source with 1 on each site for pion 
creation operator, and a strange quark source with e lp x 
on each lattice site for kaon creation operator [30| ]. We 
write the K* — > ttK and ttK — > K* three-point functions 
in terms of the quark propagators G, 

C^K'(p;h,t 2 ,t 1 )= e ip - X3 (Tr[G t2 (x 3 ,t 3 ) 7 5 

X3.X2 

x G f ti (x 3 , t 3 )~f 3 G tl (x 2 , i 2 )]75), 
C K *^K(p;h,t 2 ,ti)= e 4p - X2 (Tr[G t3 (x2,t 2 )73 

X2,X 3 

xG t t i (x 2 ,t 2 )7 5 G il (x 3 ,t 3 )]75),(19) 



5 Private communication, C. DeTar (2012). 



where trace is over the color index; the Dirac matrix are 
used as an interpolating field for the ith meson: 75 for 
pseudoscalars and 73 for the vector meson. 
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D. Extraction of energies 



B. Computations 



To map out "avoided level crossings" between the K* 
resonance and its decay products, it is important to use 
the variational method [36[ to separate the ground state 
from the first excited state by calculating a 2 x 2 cor- 
relation function matrix C(t) denoted in (fT3")) . For this 
purpose, we construct a ratio of the correlation function 
matrices as 



M(t,t R )=C(t) C-\t R ), 



(20) 



with some reference time t R [36J to extract two lowest 
energy eigenvalues E n (n = 1,2), which can be obtained 
by a cosh-fit to two eigenvalues X n (t,t R ) (n = 1,2) of 
the correlation matrix M(t,t R ). Considering the use of 
the staggered fcrmion, it is easy to verify that X n (t,t R ) 
(n = 1,2) explicitly has an oscillating term (oT| - [69| . 



X n (t, t R ) 



An cosh 



-{-lfB n cosh 



t 



T 
2 

-E' 1 1 



T 



(21) 



for a large t, which mean « tj < i C T/2 to sup- 
press both the excited states and wrap-around contribu- 
tions [30l. l7"oT - [72| 6 . Without loss of generality, we suppose 
Xi(t,t R ) > X 2 (t,t R ). 



III. LATTICE CALCULATION 
A. Simulation parameters 

We use the MILC full QCD gauge configurations in 
the presence of the Nf = 2 + 1 flavors of the Asqtad- 
improvcd staggered fermions [60l Ir3ll | and a Symanzik- 
improved gluon action (73l |. We should keep in memory 
that the MILC gauge configurations are generated using 
the staggered formulation of lattice fermions [74| with 
the fourth root of the fermion determinant (60| . 

We measured the ttK four-point correlation functions 
on the 0.15 fm MILC "medium" coarse lattice ensem- 
ble of 400 20 3 x 48 gauge configurations with the bare 
quark masses am u d = 0.00484 and am s = 0.0484 and 
bare gauge coupling 10/ g 2 = 6.566. The inverse lat- 
tice spacing a -1 = 1.373^4 GeV and the lattice extent 
L ps 3.0 fm [H, [33[. The mass of the dynamical strange 
quark is quite close to its physical value, and the masses 
of the u and d quarks are degenerate [32L l33j . The more 
detailed descriptions of the simulation parameters can 
be found in Refs. HH, [33[. The PBC is imposed to three 
spatial directions and temporal direction. 



We employ the standard conjugate gradient method to 
achieve the necessary matrix element of the inverse Dirac 
fcrmion matrix to compute the ttK four-point functions. 
We compute the correlators on all the time slices, and 
explicitly combine the results from each of the Nt = 48 
time slices. To be specific, the diagonal correlator Cn(t) 
is measured through 

C n (i) = (frK)(t)(*K)H0)) 

After averaging the propagator over all Nt = 48 possi- 
ble values, we found that the statistics are significantly 
improved. 

For each time slice, six Dirac fermion matrix inversions 
are needed to compute for the possible 3 color choices 
for the pion source and kaon source, respectively. So, to- 
tally we carry out 288 matrix inversions on a single gauge 
configuration. This big number of the matrix inversions, 
conducted on 400 MILC gauge configurations, furnishes 
the gigantic statistics required to precisely compute the 
ttK four-point correlation functions. 

For the another diagonal correlator 6^22^), K* (892) 
correlator, we simply exploit the available propagators 
measured in our previous study [28| to calculate the 
if* (892) correlator 



'(ts 



where we sum the correlator over all the time slices and 
average it too. 

We evaluate the first off-diagonal correlator C21 (t): the 
txK — > K* three-point function, through 

C 21 (t) = (K*(t)(nK)H0)) 

where the summation is over all the time slice. Through 
the relation Ci 2 (t) = CJiW) we can gratuitously gain 
the second off-diagonal correlator Ci 2 (t): the K* — > ttK 
three-point function. 

In the present study, we evaluate two-point correlation 
functions for pion and kaon as well, 

G„(0;t) = i5^<0|7rt(0,t + i,MO,t 8 )|0>, 
G K (p;t) = ^J2(0\K^p,t + t s )K(p,t s )\0), (22) 



In Ref . [30| , we gave a detailed discussion about a contamination 
from "wraparound" effects. In practice, we will select the fitting 
time ranges satisfying t m ax < 16, and reasonably neglect it. 



where the summation is over all the time slice, and the 
Gir(0; t), Gk(P', t) are the two-point correlation functions 
for pion meson with zero momentum, and kaon meson 
with the momentum p, respectively 
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IV. SIMULATION RESULTS 

In our previous work [28[ , we have measured the point- 
to-point pion and kaon correlators with high accuracy. 
Exploiting these correlators, we can precisely derive the 
pion mass (m v ) and kaon mass (thk), which are in fair 
agreement with the previous MILC determinations in 
Ref. [Hj]. In Table U we list the pion mass m v , the mass 
rriK and energy Ek of kaon meson with the momentum 
p = (27r/L)e3, which are extracted through a single ex- 
ponential fit ansatz to G v (t; 0) and G/f (i; p) in Eq. (|22[). 
We show the mass rriK* and energy Ek* of the vector 
K* meson with the momentum p = (27r/L)e3 as well, 
which are extracted from the K* (892) correlator. 



TABLE I: Masses m of the pion, kaon and A* (892) mesons, 
and energies E of kaon and A* (892) mesons with the momen- 
tum p = (27r/L)e3, extracted from the corresponding point- 
to-point correlation functions. 





IT 


K 


A* (892) 


am 


0.17503(17) 


0.39913(27) 


0.7757(70) 


aE 




0.50465(48) 


0.8278(82) 



We must stress at this point that, in the present work, 
we just use this calculated A* (892) mass rriK* to indicate 
the position of free A* (892) mass, which are marked by 
the fancy cyan plus point in Fig. and we visualize this 
value to compare with the resonance mass Mr. 



A. Diagrams D, C, and R 

The individual ratios R x (X = D,C and R), which 
correspond to the diagrams in Fig. [T] arc illustrated in 
Fig. [3]as the functions of time separation t. The values of 
the direct amplitude ratio R D are pretty close to unity, 
implying a quite slight interaction in this channel. On the 
other hand, the crossed amplitude ratio R c increases lin- 
early, hinting a repulsive force in this channel. Moreover, 
after a starting increase up to t ~ 4, the rectangular am- 
plitude ratio R R demonstrates a roughly linear decrease 
up until £ ~ 15, and the signals become noisy after that, 
suggesting an attractive force between pion and kaon in 
this channel. These characteristics are what we expected 
from the theoretical predictions [39|, [75[ . 

We can observe that the crossed and rectangular ampli- 
tudes take the same value at t = 0, and the similar values 
for small t. Since our analytical representations for both 
amplitudes are identical at this value of t, they should 
manifest analogously until the asymptotic 7rA state is 
reached. Clear signals observed up to t = 15 for the 
rectangular amplitude demonstrate that the technique of 
the moving wall source without gauge fixing used here is 
practical and feasible. 

According to the analytical arguments in Ref. (7o| , we 
can readily infer that the ratio for the rectangular di- 



0.30 
0.25 



0.20<E^><&<J > <!> < S><J > O<J><I>OOO<I> < I > O < I >< l > ^ > # < §"5'^^^ 
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0.05 

-o.oo[§QC2 
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R c (t) 



OGOOO 



R R (t) 



[]■ 



[]_ 



[] 
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t 

FIG. 3: (color online). Individual amplitude ratios R x (t) 
(X — D,C and R) for the irK four-point correlation func- 
tion measured by the moving wall source without gauge fix- 
ing as the functions oft. Direct diagram (magenta diamonds) 
shifted by 0.8, crossed diagram (red octagons) and rectangu- 
lar diagrams (blue squares). 



agram R R has errors, which should increase exponen- 
tially as e™ K * for large time separation. The magni- 
tude of the errors is in quantitatively agreement with 
this theoretical prediction as displayed in Fig. [4] Fit- 
ting the errors SR R (t) by a single exponential fit ansatz 
SR R (t) ~ exp(/z fi i) over the range 10 < t < 16, we 
can achieve the corresponding fitting values of fj,R with 
afiR = 0.358, which can be reasonably compared with 
the corresponding kaon masses ttlk determined in our 
previous work [281 ] , which is also listed in Table |TJ This 
demonstrates, on the other side, that the technique of 
the moving wall source without gauge fixing used in this 
work is practically feasible. 




FIG. 4: (color online). The error of ratio R R (t) as a function 
of time separation t. Solid line is a single exponential fit over 
the range 10 < t < 16. 
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B. Energy eigenvalues 

We calculate two eigenvalues X n (t, tn) (n = 1,2) for 
the matrix M{t,tg) denoted in Eq. (|2H)) with the ref- 
erence time tji = 5. In Fig. [5] we illustrate our lattice 
simulation results for X n (t,tfj)(n = 1,2) in a logarith- 
mic scale as a function of time separation t along with a 
correlated fit to the asymptotic form offered in Eq. (f2~T|) . 
From these fits the desired energies E n (n = 1,2) are 
then obtained, and which will be employed to derive the 
p-wave scattering phase. 




2 4 6 8 10 12 14 16 18 20 22 24 
t 



U. t ( 


d 1 


i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i i i 
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FIG. 6: The effective energy plots, aE n (n = 1,2), as the 
functions of t m i n . (a) The effective energy plot for E\ and (b) 
that for E'z 



FIG. 5: (color online). The eigenvalues Ai(i, tn) and \2{t, tn) 
as a function of t. Occasional points with negative central 
values for the eigenvalue \2{t,tn) are not plotted. The solid 
lines are the correlated fits to the asymptotic form denoted in 
Eq. (|21|l. from which the energy eigenvalues E n (n = 1, 2) are 
extracted. The lower curve (n = 2) is slightly steeper than 
the upper curve (n = 1). 

As we noticed in Refs. [29l - [3l| . we realize that the prop- 
erly extracting the energy eigenvalues is vital to our final 
conclusions. Since the PBC is imposed on three spatial 
directions and the temporal direction, we should suppress 
the warp-around contaminations [Til l3l| . By defining a 
fitting range [i m i n ,£ max ] and varying the values of the 
minimum fitting distance t m i n and the maximum fitting 
distance i max , we obtain these energies in a correct man- 
ner. In practice, we make t m j n = tn + 1 and increase ref- 
erence time tfi to reduce excited contaminations fll| . At 
the same time, we opt £ max to be away from the time slice 
T/2 to reduce the warp-around effects [ll|. Furthermore, 
we extract two eigenvalues A n (n = 1,2) with the "effec- 
tive energy" plots [ll], Ull > a variant of the effective mass 
plots, and they were fit to Eq. ([21]) by changing t m i n , and 
with the t max either at 15 or where the fractional statis- 
tical errors exceeded about 20% for two successive time 
slices. The effective energy plots as a function of i m ; n arc 
illustrated in Fig. [5] 

The energy eigenvalues E„ (n = 1,2) were chosen by 
looking for the combination of a "plateau" in the effective 



energy plots as the function of t min and a reasonable fit 
quality [111, [3l| . We observed that the effective energies 
take only relatively small errors within a minimum time 
distance region 5 < t min < 8 for E\ and 5 < t min < 6 for 
E2, respectively. The fitted numbers for E n (n = 1,2) 
along with the fitting parameters tn, t m - ln and £ max , fit 
quality x 2 /dof are summarized in Table Ull 

TABLE II: The fitted values of the energy eigenvalues for the 
ground state (n = 1) and the first excited state (n — 2). 
Here we tabulate the reference time tn, the fitting range, t m i n 
and t max , the number of degrees of freedom (dof) for the fit 
quality x 2 /dof and the fitted results for the energy eigenvalues 
E n (n — 1, 2) in lattice units. 



n 


t R 


^min 


^max 


aE n 


X 2 /dof 


1 


5 


6 


15 


0.67507(40) 


12.2/6 


2 


5 


6 


15 


0.8534(78) 


9.6/6 



The energy of the pion and kaon in non- interacting case 
(namely, E\) is computed from the pion mass and 
kaon energy Ek listed in Table Q] as E\ = iu^ + Ek- This 
number is listed in Table InTI We distinctly observe that 
E\ < E\ < E2, which mean that the p-wave scattering 
phase for Xi(t,tu) and A2 (t,tjj) is positive and negative, 
respectively. This evidently reveals that there exist a 
resonance between E\ and £"2- 
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TABLE III: Summary of the energy eigenvalues E n (n = 1, 2) 
and p-wave scattering phase Si for irK system in a cubic box. 
Ei is the energy of the free pion-kaon system. E n (n = 1,2) 
are obtained from the correlated fits to eigenvalues \ n (t,tn) 
(n = 1, 2). The invariant mass y/s, the scattering momentum 
k and the p-wave scattering phase Si derived through the 
expression ([23} in the continuum are regarded as Cont, and 
those achieved with the equation (|24p on the lattice are viewed 
as Lat. The scattering momentum ko is computed by kg = 
1/4 x (y/s + (m 2 — rn 2 K )/yfs) 2 — m 2 . All values with the mass 
dimension are in lattice units. 





n=l 




n = 2 




E n 


0.67968(51) 








E n 


0.67507(40) 




0.8534(78) 






Cont 


Lat 


Cont 


Lat 




0.59751(45) 


0.60350(45) 


0.7934(84) 


0.8004(84) 


k 2 


0.00588(13) 


0.00750(14) 


0.0690(30) 


0.0729(31) 


h 2 




0.00745(13) 




0.0717(30) 


tan<5i 


0.0294(93) 


0.0091(22) 


-2.01(43) 


-2.48(63) 


sin 2 Si 


0.00087(55) 


0.000083(39) 


0.802(68) 


0.860(62) 




C. 


Finite-size 


effects 





Wc should pay attention to the discretisation error (or 
truncation error) inherent in the special finite size for- 
mula for nK system denoted in (fTTj) [2lL |30| . It stems 
from the Lorentz transformation from the moving frame 
to the center-of-mass frame. When applying the Lorentz 
symmetry in the continuum limit, we utilize the following 
relations [HIH, 



k z = 



1 K 



(23) 



in the Lorentz transformation for the invariant mass yfs, 
the energy of the ttK system in the moving frame Emf 
and the scattering momentum k. Nevertheless, the dis- 
cretisation effects definitely violate the Lorentz symme- 
try on the lattice and equation (|23j) is only effective up 
to the truncation errors. 

Following the recommendations in Ref. [3(|, we calcu- 
late the invariant mass yfs and the scattering momentum 
k from the energy in the moving frame Emf of nK sys- 
tem using 



cosh(\/i) = cos!i(£'mf) - 2 sin 2 (J^j 



2sin 2 (fc/2) = cosh I ^~ 



' K 



2yf S 



- cosh(m T ), (24) 



and derive the p-wave scattering phase 5\ by inserting 
the scattering momentum k into the finite-size formula 
in Eq. (fTTj) . We have justified these formula in Ref. [30l ]. 
and we will employ them in this work to investigate the 
discretisation effects. 



To grasp these discretisation effects quantitatively, in 
the present study we compute the invariant mass yfs and 
the scattering momentum k not only from the energy mo- 
mentum relation in the continuum (|23p but also from that 
on the lattice (f24j) . and subsequently extract the p-wave 
scattering phase S\. People usually regard the dispar- 
ity stemming from two options of the energy momentum 
relations as the discretisation error. Since it is a kind 
of truncation error, it is expected to smaller if we em- 
ploy the lattice ensemble with smaller lattice space a, of 
course it should be disappeared in the continuum limit. 
The results for the p-wave scattering phase 5\ along with 
the invariant mass yfs and the scattering momentum k 
are summarized in Table [ 



D. 



Extraction of the scattering phase and decay 
width 



The noticeable differences in the invariant mass yfs 
and scattering momentum k because of the discretisation 
effects are obviously observed from Table Mil Moreover, 
the differences for the p-wave scattering phase 5\ due 
to the discretization effects are impressive, and can be 
comparable with the statistical errors, even considerably 
larger than its statistical error for the n = 1 case. These 
characteristics are visualized in Fig. where the p-wave 
scattering phase sin 2 Si is displayed instead [HI, |s3|. We 
notice, from Table \LU\ that the numerical value of the p- 
wave scattering phase Si at the invariant mass yfs < rriK* 
(arriK' = 0.7757±0.0070) is positive due to an attractive 
interaction, and that at yfs > uik* is negative owing to a 
repulsive interaction. These features indicate a resonance 
at a mass around the K* mass. 

In principle, it is a piece of cake to extract the K* (892) 
meson decay width through fitting the p-wave scatter- 
ing phase shift data with the effective range formula di- 
rectly [l6|,[3(|. Nevertheless, in this work we studied with 
the quark mass which is larger than its nature value and 
the kinematic factor in the decay width depends clearly 
upon the quark mass p7| . thus an extrapolation is in- 
dispensable [HI, [3(|. However, we carried out a lattice 
calculation with one set of the quark mass in this ex- 
ploratory investigation, therefore, we have no choice but 
to adopt an alternative method [3, H(| . As we explained 
in section III Al the resonant characteristic of the p-wave 
scattering phase S\ are parameterized with the coupling 
constant gn'-nK, 



tan^i = 



9k*tvK 



6ir y/s(Mft 



(25) 



where Mr is the resonance mass. 

According to the elaborations in Refs. [13, [z3, we can 
fairly suppose that the coupling constant gx'TrK changes 
quite slowly and smoothly with the quark mass. There- 
fore, the equation (j2"5j) enables us to solve for two un- 
known quantities, that is, the coupling constant gK*irK, 
and the resonance mass Mr [lj| [3(J • 
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According to the discussions in Rcfs. [Ty, [30|, in prac- 
tice, we usually employ the scattering momentum fc 
instead of k when applying Eq. (|25|) . In Table IIII| 
we provide the scattering momentum fco calculated by 
fco = 1/4 x + (?n 2 — m 2 K ) / yfs) 2 — m 2 in addition to 
k. We can observe that the difference between k and ko 
is not significant, and we can neglect this systemic error 
for the present study [HI, H(| ■ 

When we utilize the energy-momentum relations (f23|) 
in the continuum, the lattice simulation results of the 
coupling constant qk'-kK and the resonance mass Mr 
solved by Eq. {25} are 



9K"nK 

Mr 
M R /m K * 



11.73 ±2.08, 

0.739(20), 

0.953(28), 



(26) 



where the K* meson mass mx* is quoted from our pre- 
vious work [28j . 

On the other hand, if we adopt the energy momen- 
tum relations (|24l) on the lattice, we gain the simulation 
results as 



g K **K = 6.38(78), 
M R = 0.7873(97), 
M R /m K , = 1.015(16). 



(27) 



This obtained value of the coupling constant qk*-kK in 
lattice case is in fair agreement with gK*nK ~ 5.5, which 
are obtained by Nebreda and Pelaez from the residue of 
the amplitude at the pole position in Ref. (27j. More- 
over, it is in reasonable agreement with the experimental 
observable gn'-nK = 5.64(35) evaluated from the PDG 
estimations of the decay width IV* = 50.8(9) MeV [J 
within the statistical error. 

In Fig. [3 we illustrate the curves for sin 2 <5i achieved by 
equation (|25[) with the coupling constant qk-^k and the 
resonance mass Mr provided in Eq. (f2"6"|) and Eq. (|2"T|) . 
respectively. The positions of the resonance mass Mr, 
are courteously marked in Fig. [7] for the two cases (black 
cross and red plus for the continuum and lattice cases, 
respectively). For visualized comparisons, we mark the 
K* (892) mass ttik* with fancy cyan plus as well. We can 
observe that the resonance mass Mr for lattice case is in 
reasonable accordance with the if* (892) mass xtik*- 

Supposing that the dependence of the coupling con- 
stant gK*7rK on the quark mass is ignorable [27], [77], we 
can roughly estimate the K* (892) meson decay width at 
the physical quark mass as 



r phy _ 9 2 K^K (fc PhY ) 



6^ [mfiy ' 



(28) 



where rn^I = 891.66(26) MeV is the physical if*(892) 
meson mass, which we take from the current PDG [l|, 
and the scattering momentum fc phy at the physical point 
is calculated by 



( fcP hy)5 



1 



,phy 



phy\2 ' 



(mf0M_mf) 
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FIG. 7: (color online). The scattering phase sin 2 Si and the 
positions of the K* meson mass itik* and resonance mass Mr. 
The simulation results achieved with the energy-momentum 
expressions in the continuum (|23[) are viewed as Cont and 
those with the relations on the lattice (|24p as Lat. The two 
lines are achieved by Eq. (|25|l with the quantities gK*-wi< and 
M R provided in Eq. and Eq. (|2"7). respectively. The 

abscissa is in lattice units. 



where mP hy is physical pion mass (mP hy = 139.57018(35) 



MeV) and m^ y is physical kaon 



mass 



493.677(13) MeV) which are quoted from the PDG 0]. 
This produces 



pphy = ( 219 ± 39) MeV 



(29) 



where we utilize the simulation result given in Eq. (|26l) . 
On the other hand, if we use these given in Eq. (|2"T|) . it 

yields 



r phy = ( g4 9 ± 8 ) Me y_ 



(30) 



The estimation in lattice case given in Eq. (|3"D"|) is in fair- 
agreement with the corresponding PDG quantity for the 
K* ttK decay width, namely, T K » = 50.8 ± 0.9 MeV. 
We can observe that the difference stemming from our 
two options of the energy-momentum relations is much 
larger than the statistical error, which indicate that the 
proper correction of the discretization errors is absolutely 
necessary. 

These results are quite stimulating, considering that we 
make a big assumption that the coupling constant gx*-KK 
is independent upon the quark mass, and we carry out an 
extrapolation, etc. Anyway, one thing greatly comforts 
us is that we use the pion mass (about 240 MeV) which 
is pretty close to its realistic value (about 140 MeV), so 
we don't carry out a long extrapolation. 
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V. CONCLUSIONS AND OUTLOOKS 

In the present work, we have carried out a lattice QCD 
computation of the p-wave irK scattering phase in the 
7 = 1/2 channel near the if* (892) resonance region in 
the moving frame with total non-zero momentum for the 
MILC "medium" coarse (a ps 0.15 fin) lattice ensemble 
with the Nj = 2 + 1 flavors of the Asqtad improved 
staggered fermions. We employed the technique with the 
moving wall source without g aug e fixing (3l| introduced 
by Kuramashi et al. in Rcfs. 140. |4 ll to calculate all the 
three diagrams classified in Ref. |43j with high precision. 
We have exhibited that the lattice QCD computation of 
the p-wave scattering phase for the 7=1/2 irK system 
and then the estimation of the decay width of if* (892) 
meson are feasible even with our present limited comput- 
ing resources. The scattering phase data clearly reveals 
a resonance at a mass around the if* (892) meson mass 
obtained in our previous study (28[. Most of all, we ex- 
tracted the if* (892) meson decay width from the scat- 
tering phase data and demonstrated that it is reasonably 
comparable with the if* (892) meson decay width quoted 
from PDG within the statistical error. 

We have adopted the effective range formula, which 
allows us to exploit the effective if* — > irK coupling con- 
stant gK'irK to extrapolate from our lattice simulation 
point (m w + TOif)/mjf* = 0.7388 to the physical point 
(tott + 171k) /ma* = 0.7102, assuming that the coupling 
constant gK*nK is independent of the quark mass. This 
is a crude estimation, a more rigorous computation of the 
decay width is highly desirable. As we pointed out above, 
the decay width can be reliably estimated from the en- 
ergy dependence of the scattering phase data by fitting 
the BWRF if we carry out the lattice simulations close to 
the physical quark mass and obtain several energies near 
the resonance mass. We will keep on enthusiastically 
requesting for the possible computational allocations to 
fulfil this valuable work. 

Nevertheless, we should bear firmly in mind that some 
critical issues should be resolved in the more convincing 
calculation. One is to reduce the discretization errors, 
which, we illustrated in the previous section, are signif- 
icantly larger than the corresponding statistical errors. 
A naive way to handle this question is to utilize a lat- 
tice ensemble closer to the continuum limit. Another 
challenging and stimulating topic is to suppress the con- 
taminations of the p-wave scattering phase from the d- 
wavc scattering phase or higher, which we prelimin arily 



touched on this topic for the 7rif system in Ref. [21 
see more valuable discussions in Ref. [22j. Moreover, the 
comprehensive investigations on the lattice size depen- 
dence of the scattering phase is absolutely fascinating. 
Nevertheless, all of these open questions are beyond the 
scope of this paper since this will demand a huge amount 
of computing allocations. We postpone these expensive 
tasks in our future study. 

This work concentrated on the p-wave scattering phase 
only at two energies for a single lattice ensemble. Since 



we had only a small number of energies at hand, it be- 
comes quite difficult to reliably map out the resonance 
region. Therefore, our current lattice results are not 
comparable with the experimentally measured quantities. 
Although a reliable derivation of the if* (892) resonance 
parameters from the lattice is absolutely big challenging 
and most prospective, our rudimentary work reported 
here can be still viewed as an important conceptual study, 
and the techniques employed here will be helpful and use- 
ful for other resonances such as the D* , possibly even for 
some exotic hadrons. 
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Appendix A: The numerical evaluation of the 
Z2o(l;q 2 ) function 

In this appendix we follow the original derivations and 
notations in Refs. |22l |30L l37| to provide one simple ap- 
proach for the numerical evaluation of the zeta function 
Z$ Q (s; q 2 ) defined in Eq. JUJ in the moving frame for the 
arbitrary value of q 2 . 

The definition of the zeta function Z^s; q 2 ) appeared 
in Eq. flU]) is 



Z^a 2 ) 



(Al) 



where 3^m(r) = r l Yi m (Sl r ). f2 r represents the solid angles 
(8, <j>) of r in spherical coordinates and the Yi m are the 
spherical harmonic functions, and the summation for r is 
taken over the set 



= {i 



r = 7 



n 



where 



a = 1 + 



2 



n e Z 3 } , (A2) 



E 2 

11 'CM 



and the operation 7 1 is defined in Eq. (J5J. 



13 



Without loss of generality, we first assume q 2 > 0, and where the summation over r is conducted with rgPj de- 
divide the summation in Z 2 q into two pieces as noted in Eq. (|A2|) . The second term can be conveniently 



E 



3^20 (r) 3^20 (r) v-^ 3^20 (r) 



( r 2_ ? 2) S ( r 2 2) S '^ ( r 2_ g 2) S 

rGP d V H ' r 2 <q 2 V H ' r 2 >q 2 V y ' 



(A3) 



delivered in an integral expression, 



3^20 (r) 



r 2 >q 2 y i 



ch t-i e -*^-«") 

3^20 (r) 



rEPd 



(r 2 -q 2 ) 1 



E ^ (r) J 



(r 2 - g 2 ) s 



r 



(A4) 



first term 



The second term nicely counteract the first term in 
Eq. (|A3|) . Using the Poisson resummation formula, the 
first term results in 

rh) f dtts ~ 1&tq2 E 

/„ = J d 3 x3; 2 „(r)e- t l r l 2+ ' ;2 ™- x , (A5) 

where r = 7 _1 (x+iad). After transforming the integra- 
tion variable from x to r, and considering the relations: 
rf 3 x = r yd 3 r and x = ■yr — ^ad, then we can separate 
terms which are dependent only upon r 

j n _ 7e -i7r«n.d J d 3 r -y 2 ^ r - )e -«|r| 2 +i27r7n.r_ 

Let k = ir-yn, we rewrite above equation as 



where r = (x,y,z). Let us conduct a variable substitu- 
tion, namely, r — ik/t — > r, we can strictly verify 



d 3 rJ/ 2 e -t(r-ik/f) 

We finally obtain 



2k 

T 
7 

7 

7T 

7 



dx y 


x 2 - 


(h 


kx 




t 2 




k 2 




t 2 


(h 


k 2 




t 2 



h 2s 

h x \ -tx z 

t 2 



(A6) 



/„ = _ 7 e-" rQn ' d e~ fc2/ * 2 ' 



3/2 



t 7/2 



3>2 (k), 



where 3^20 (k) = fc 2 ^2o(^fc)- Now we can rewrite the first 
term in Eq. (|A4[) as 



first term 



1 



r(s) Jo «* 



xe i7ran-d e -(i7T7n) /t 



(A7) 



After collecting all terms, it leads to the representation 
of the zcta function at s = 1, 



^ d o(i;? 2 ) 



£ r 2 F 20 (r! r )' 



r£P d 



r 2 — q 2 



die 1 " — V (7r 7 n) 2 r 20 (^ 



xe -iiran-d e -(7r7n) /t^ 



(A8) 



Now let us consider the case of q 2 < 0. As a matter 
of fact, there is no need for us to divide the summation 
in Zf (s;q 2 ) into two parts, and can be delivered in an 
integral expression as well [30l . l3lj . Conducting the same 
procedures, we obtain the same expression in Eq. (|A8[) . 
Therefore, equation (|A8[) is applicable for the arbitrary 
value of q 2 . 

Plugging in d = (0, 0, 1) into equation (|A8[) . we arrive 
at the representation of the zeta function Zf (s;q 2 ) in 
Eq. (TT2]) devoted for the current work 



2 2 d (i;g 2 ) 



= r2y 2o(0-)' 



r£P d 



r 2 — q 2 



die 4 ? — V (7T7n) 2 r 20 (^ 
'° l ~ 2 neZ 3 

xcos(7ran-d)e-^ n ) 2 /*, 



(A9) 



where only the real part of the zeta function Z d (l; g 2 ) 
is survived. 

I also note that the general numerical evaluation of the 
zeta function Zf m (s\ q 2 ) has been derived in Refs. [2lll22l ]. 
We numerically compared both these representations of 
the zeta function iJ 2 o(l; q 2 ) with this representation, and 
found that they are reasonable consistent. 
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